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We study the big-bang nucleosynthesis (BBN) scenario with late-decaying exotic particles with 
lifetime longer than ~ 1 sec. With a late-decaying particle in the early universe, predictions of 
the standard BBN scenario can be significantly altered. Therefore, we derive constraints on its 
primordial abundance. We pay particular attention to hadronic decay modes of such particles. We 
see that the non-thermal production process of D, 3 He and 6 Li provides a stringent upper bound on 
the primordial abundance of late-decaying particles with hadronic branching ratio. 
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It has long been recognized that the (standard) big- 
bang nucleosynthesis (BBN) provide a good probe for 
the early universe. With our current knowledges of nu- 
clear reaction processes, we can precisely calculate abun- 
dances of the light elements (in particular, D, 3 He, 4 He, 
6 Li, and 7 Li) as functions of the baryon-to-photon ra- 
tio r\ = ub/tIj. Thus, comparing the theoretical pre- 
dictions with the observations, we can obtain various 
informations about the evolution of the universe. Im- 
portantly, using the value of rj suggested by the WMAP 
(t) = (6.1 ± 0.3) x 10~ 10 [1]), the theoretical predictions 
show relatively good agreement with the observations. 

In cosmological scenarios in the frameworks of physics 
beyond the standard model, however, the BBN may not 
proceed in the standard way. This is because, if we as- 
sume physics beyond the standard model, there exist var- 
ious exotic particles. Those exotic particles may cause 
non-standard processes and spoil the success of the stan- 
dard BBN. 

In particular, if the exotic particle (called X hereafter) 
decays radiatively and/or hadronically after the BBN 
starts, the primordial abundances of the light elements 
may be significantly affected. Indeed, energetic particles 
produced by the decay of X may scatter off and dissoci- 
ate the background nuclei. If such processes occur with 
sizable rates, predictions of the standard BBN scenario 
are changed. 

In various models of particle physics, there exist long- 
lived (but unstable) particles and hence their effects on 
the BBN should be studied. In particular, many of those 
particles interact very weakly and hence it is difficult to 
study their properties by collider experiments. Thus, in 
some case, the BBN provides useful and important infor- 
mations about those weakly interacting particles. Proba- 
bly the most famous example of such late-decaying parti- 
cle is the gravitino in the supergravity theory. Gravitinos 
are produced in the very early universe by scattering pro- 
cesses of the particles in the thermal bath. Their interac- 
tion is suppressed by the inverse powers of the (reduced) 
Planck scale M», and hence the lifetime becomes very 
long. An order of magnitude estimate shows that, if the 
gravitino is lighter than ~ 10 TeV, its lifetime becomes 
longer than ~ 1 sec. In this case, the thermally produced 



gravitinos decay after the BBN starts. 

Effects of the radiative decay of such long-lived parti- 
cles have been extensively studied (see, e.g., [2-5]). How- 
ever, in many cases, hadronic branching ratio may not 
be negligible and hence, in studying the BBN with late- 
decaying particles, it is necessary to consider effects of the 
hadronic decay processes. Even if X dominantly decays 
into the photon (and something else) , hadronic branching 
ratio is expected to be at least 10 - ( 2-3 ) since the emitted 
photon can be converted to a quark- antiquark pair. Of 
course, if X directly interacts with the colored particles, 
hadronic branching ratio may become larger. 

In the past, the BBN with hadro-dissociation processes 
induced by hadronic decays of long-lived particles was 
studied in [6], which are effective for relatively long life- 
time (> 10 2 sec).^ 1 The analysis in Ref. [6] contains, 
however, a lot of room to be improved since many of nu- 
clear reactions that they used were not accurate enough 
or not available at that time. After the study of [6], 
however, there have been significant theoretical, experi- 
mental and observational progresses in the study of the 
BBN. First of all, new data for the hadron reactions have 
become available and their qualities have been improved 
very much. Moreover, the primordial abundances of the 
light elements have been precisely determined with var- 
ious new observations. In addition, it has been recently 
known that some of the non-standard processes induced 
by the decay of late-decaying particle X, which were not 
taken into account in [6], may play important roles in 
the BBN. With these progresses, a new study of the late- 
decaying particles with hadronic branching ratio should 
be relevant. 

Thus, in this letter, we reconsider the BBN processes 
with long-lived exotic particle X paying particular at- 
tention to the effects of the hadronic decay modes. As a 
result, we will see that, with hadronic decay modes, the 
constraint on the primordial abundance of X becomes 



#1 BBN constraints from the interconversion process between 
neutrons and protons by hadronic decays were studied in 
Refs. [7,8], which is effective for shorter lifetime (< 10 2 sec). 
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very severe compared to the case only with the radiative 
decay modes. In particular, we will see that non-thermal 
production of D, 3 He and 6 Li provides a stringent con- 
straint. 

We first introduce the framework of our study. 
Although we have several candidates of the late- 
decaying particles, we perform our analysis as model- 
independently as possible. Thus we parameterize the 
property of the late-decaying particle X using the fol- 
lowing parameters: ex (released energy from the single 
decay of X) which is equal to its mass mx unless specially 
stated, Ejet (energy of the primary parton from the de- 
cay of A), tx (lifetime), (hadronic branching ratio), 
and the primordial abundance of X. We parameterize 
the primordial abundance by using the following "yield 
variable" Yx = nx / s, which is defined at the cosmic time 
t -C Tx- Here, nx is the number density of X while s is 
the total entropy density. We assume that X decays only 
into the particle in the "observable sector," and that the 
branching ratio decaying into the hidden-sector particle 
vanishes. 

In our analysis, we calculate the primordial abun- 
dances of the light elements for given sets of the parame- 
ters listed above taking account of dissociation processes 
due to hadronic (and electromagnetic) interactions. The 
outline of our calculation is as follows. In order to study 
the effects of the hadronic decay modes, we first calculate 
the energy spectrum of the (primary) hadrons (in partic- 
ular, protons and neutrons) generated from the partons 
directly emitted from X. Then, we follow evolutions of 
the hadronic showers. As a result, the numbers of the 
light elements produced (or destroyed) by the decay of 
X is calculated. We include the hadro-dissociation pro- 
cesses (as well as the photo-dissociation ones) into the 
BBN reaction network and numerically follow the evolu- 
tion of the abundances of the light elements. Details of 
our study will be described elsewhere [9]. 

When X decays hadronically, quarks or gluons (i.e., 
partons) are first emitted. Those partons are hadronized 
soon after the decay. Thus, we should consider how 
the nucleons (and the mesons) propagate in the thermal 
bath. 

In considering the propagation of stable particles (i.e., 
proton, neutron, and heavier nuclei), there are two classes 
of important processes. The first is the scattering process 
with the background photons 7bgS and electrons eg G s. 
By scattering off 7bg s an d eg G s, energetic nuclei lose 
their energy without affecting the abundances of light 
elements. The second class is the scatterings with the 
background nuclei. With such processes, first of all, the 
background nuclei become energetic after the scattering. 
Therefore the energetic nuclei are copiously produced. 
(We call this "hadronic shower.") In addition, if inelas- 
tic scattering occurs, background nuclei are dissociated 
and the abundances of the light elements are changed. 
Thus, if the second class of reactions occur significantly, 
the abundances of the light elements deviate from the 



predictions of the standard BBN. 

Including relevant hadronic scattering processes (as 
well as photo-dissociation processes), we have calculated 
the abundances of the light elements. In our study of the 
evolution of the hadronic showers, the basic framework 
is the same as that used in [6] although there are several 
modifications. The most important improvements are as 
follows, (i) We carefully take into account the energy 
loss processes for high-energy nuclei through the scatter- 
ing with background photons or electrons. In particular, 
dependence on the cosmic temperature, the initial ener- 
gies of nuclei, and the background 4 He abundance are 
considered, (ii) We adopt all the available data of cross 
sections and transfered energies of elastic and inelastic 
hadron-hadron scattering processes, (iii) The time evo- 
lution of the energy distribution functions of high-energy 
nuclei are computed with proper energy resolution, (iv) 
The JETSET 7.4 Monte Carlo event generator [10] is 
used to obtain the initial spectrum of hadrons produced 
by the decay of X. (v) The most resent data of observa- 
tional light element abundances are adopted, (vi) We es- 
timate uncertainties with Monte Carlo simulation which 
includes the experimental errors of the cross sections and 
transfered energies, and uncertainty of the baryon to pho- 
ton ratio n. (We take n = (6.1 ± 0.3) x 10~ 10 .) 

We are interested in the situation where the number 
density of X is small enough so that the energetic par- 
ticles in the hadronic shower dominantly scatter off the 
background particles. (Otherwise, the abundances of the 
light elements are so affected that the results become 
inconsistent with the observations.) If the energetic nu- 
clei lose most of their energy by scattering off 7bg s an d 
eg G s, the number densities of the (background) nuclei is 
almost unaffected. On the contrary, if the scattering rate 
with 7bgs and eg G s becomes negligible, number densi- 
ties of the light elements are significantly changed by the 
hadronic processes. 

For charged nuclei, energy-loss rates due to the scat- 
terings with the 7bg and eg G depend on the velocity of 
the nuclei [7] . For the temperature T < 20 keV, relativis- 
tic nuclei do not lose their energy by the scatterings with 
7bg s an d eg G s. For non-relativistic nuclei with velosity 
(3n, energy-loss is dominated by the scattering process 
with background electrons with velosity /3 e < (3n- (As 
pointed out in [7], energy loss via the scatterings with 
high- velosity electron is extremely suppressed and is neg- 
ligible in our case.) Thus, we use the energy- loss rate 

dE N ATia 2 Z 2 n e / — { Am e /3 2 N \ 

where w p is the plasma frequency, m e is the electron 
mass, Z is the charge of the nuclei, a is the fine struc- 
ture constant, A is a constant of 0(1). (In our numerical 
calculations, we take A = 1.) In addition, 

J(r) = -4= I dxx 2 e- x \ (2) 
V 71 ' Jo 
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Notice that the number density of the bac kground elec- 
tron with (3 e < (3n is given by n e I(@N / '\/2T : /m e ). For 
(3n 3> \/T/m e , I ~ 1 and the energy- loss is very ef- 
ficient. On the contrary once the velosity of the nuclei 
becomes smaller than the ther mal velo sity of the electron 
(fie) ~ y/T/me, I ~ O((0N/y/T/m e ) 3 ) and energy loss 
becomes less effective. # 2 

For the neutral particle (i.e., neutron), on the con- 
trary, the scattering process with eg G is not important 
for T < 0.1 MeV. Thus, once energetic neutrons are pro- 
duced, they may induce hadronic showers by scattering 
off the background nuclei and change their abundances. 
If the temperature becomes lower than ~ 0.3 keV, how- 
ever, most of the (energetic) neutrons decay into protons 
before scattering off the background particles. Thus, if 
Ti>3x 10 7 sec, the most important constraint on Y\ 
is from the photo-dissociation processes of the light ele- 
ments. 

In our analysis, we also take account of the photo- 
dissociation processes induced by high-energy photons 
generated from the decay products of X. Once X de- 
cays, energy as large as (1 — B h )e x is directly deposited 
into the radiation while the rest of the released energy 
Bh^x first goes into the hadronic sector. Importantly, 
however, spectrum of the high energy photon primarily 
depends on the total amount of the injected energy (as 
well as the temperature). Thus, we approximate that 
all the emitted energy in the decay process is eventually 
converted to the form of radiation, calculate the photon 
spectrum, and estimate the photo-dissociation rates of 
the light elements. For details of the treatment of the 
radiatively decaying particles, see [3]. 

Moreover, pions and kaons produced in the hadroniza- 
tion are also considered. For T > 0.1 MeV (f < 10 2 sec) 
they can scatter off the background nuclei before their de- 
cay, which leads to increase the neutron to proton ratio 
and hence the abundance of D and 4 He. In this letter, we 
have conservatively omitted the similar effects induced 
by nn and pp pairs which are effective for t > 10 2 sec [8] 
because of insufficient data for the cross sections of en- 
ergetic n and p. If we include them, the constraint can 
become severer [8] . 

With the procedure explained above, we have calcu- 
lated the abundances of the light elements for a given 
set of the model parameters. By comparing the re- 
sults with the observations, we derived upper bounds 
on Yx- As observational constraints on the light ele- 
ment abundances, we adopt the following values, D/H 



#2 In the earlier version of this letter, we did not take into 
account the velocity distribution of e^ G and used the formula 
for the case when all the background electrons have the same 
velosity (/3 e ). Thus, for /3jv < (/3 e ), the energy-loss rate we 
used was incorrect and was underestimated. This resulted in 
an overestimate of the non-thermally produced 6 Li. 



= (2.8 ± 0.4) x 10" 5 [11], 4 He mass fraction Y p = 
0.238 ± 0.002 ± 0.005 by Fields and Olive (FO) [12] and 
Y p = 0.242 ± 0.002(±0.005) syst by Izotov and Thuan 
(IT) [13], log 10 ( 7 Li/H) = -9.66 ± 0.056(±0.3) syst [14], 
6 Li/ 7 Li < 0.07(2ct) [15], and 3 He/D < 1.13(2ct) [16]. 
The above errors are at 1 a level unless otherwise stated. 
Here we added the systematic error to Yp (IT) which is 
as same as the value in Fields and Olive. For 7 Li/H we 
also added the additional systematic error for the possi- 
bilities that 7 Li in halo stars might have been depleted 
in stars or supplemented by production in cosmic-ray in- 
teractions [17,18]. For 7 Li/H we conservatively use only 
the upper bound since the experimental data for the non- 
thermal 7 Li production by energetic 4 He, especially for 
the transfered energy to 4 He through the inelastic colli- 
sions between nucleon and 4 He, are insufficient, by which 
we do not include them in our computation. 

In Figs. 1 and 2, we show the upper bound on mxYx 
as a function of the lifetime Tx for B^ = 10~ 3 and 1. 
(Here, we used mx = 1 TeV although, even if we vary 
m x within 0.1 — 100 TeV, the bounds do not change 
significantly [9].) In deriving the bound we estimate the 
confidence levels by the \ 2 fitting including both the the- 
oretical and the observational errors. As we mentioned, 
for the case with very long lifetime, upper bound is al- 
most the same as the one for the case with radiatively de- 
caying particles [3] . The most stringent constraint comes 
from the overproduction of 3 He.# 3 

For shorter lifetime t x ~ 10 2 ~ 8 sec, overproductions 
of D and 6 Li provide constraints on m x Y x . D is mainly 
produced by the hadro-dissociation of 4 He. Notice that 
the hadro-dissociation of 4 He is possible with energetic 
p and n even at such early epoch although the photo- 
dissociation of 4 He can be effective only at t > 10 6 sec 
(T < 1 keV) since high energy photons lose their en- 
ergy by scatterings off the background photons before 
interacting with 4 He at earlier epoch. The overproduc- 
tion of 6 Li is mostly due to the non-thermal processes 
with energetic T and 3 He which are generated by the 
hadro-dissociation of 4 He. These T and 3 He can be suffi- 
ciently energetic and may scatter off the background 4 He 
to produce 6 Li via the processes T + 4 He ^ 6 Li + n and 
3 Hc + 4 Hc^ 6 U+p. 

One might think that the produced 6 Li may be ener- 
getic and scatter off background nuclei to be destroyed. 
In our numerical calculations, we obtain the energy dis- 
tribution of the non-thermally produced 6 Li and calcu- 
lated the surviving rate of such 6 Li. We have found that, 
for Tx > 100 sec, the surviving rate is almost 1. There- 
fore, one does not have to take account of the destruction 
of the produced 6 Li. If the temperature is higher than 



#3 However, even for tx > 10 8 sec, the hadronic processes 
give some contribution to the D production, which can be 
seen from comparison between Fig. 1 and Fig. 2. 
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FIG. 1. Upper bounds on mxYx as a function of tx at 

95% C.L. for the case of B h = 10" 3 . The name of the 
element which gives the constraint is written by each line. 
We assume that two hadron jets are produced by single de- 
cay of X with the energy _Ej ct = mx/2. Here we consider 
mx = tx = 1 TeV. Note that Yx = nx/s. 
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FIG. 2. Same as Fig. 1, except for B b — 1. 



~ 10 keV, however, non-thermally produced 6 Li is de- 
structed by the thermal process 6 Li(p, 4 He) 3 He, as was 
also pointed out in [19]. Thus the constraint from the 
6 Li is weakened for tx ^ 10 4 sec. 

If the lifetime of X is short (i.e., Tx ^5 10 2 sec), X 
decays when the number density of the background elec- 
tron (and positron) is still abundant and, in this case, 
energetic nucleons and nuclei are likely to scatter off the 
background electron and lose their energy without com- 
mitting the hadro-dissociation processes. In this case, 
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TABLE I. Upper bounds on Tr for several values of m 3 / 2 . 

the production of D and 6 Li is suppressed and the upper 
bound on mxYx is not stringent. 

Finally, we apply the above results to the primordial 
abundance of gravitinos. In the inflationary universe, 
gravitinos are produced by the scattering processes of 
the thermal particles. The yield variable of the gravitino 
is proportional to the reheating temperature Tr after the 
inflation, Y x = 1.5 X 1CT 12 x (T R /10 10 GeV) [3]. In addi- 
tion, assuming the (massless) gauge boson and gaugino 
as the final state, lifetime of the gravitino is given by 
t 3/2 ~ 4 x 10 s sec x 7V G - 1 (m 3/2 /100 GeV)~ 3 , where N G 
is the number of the generators of the gauge group, and 
m 3/2 is the gravitino mass. As examples, we consider 
two typical cases. The one is the case where the grav- 
itino dominantly decays into the photon and photino, 
producing two hadron jets with Tj 0t = |m 3/ i 2 ; in this 
case, we take B^ = 1CP 3 , Nq = 1, and ex — h m 3/2- The 
other is the case where the gravitino dominantly decay 
into the gluon and gluino, producing one hadron jet with 
-Ejct = 5^3/2! m this case we take — 1, Nq = 8, 
and ex = ^ m 3/2 • F° r these cases, we read off the upper 
bound on the reheating temperature for several values of 
the gravitino mass. The results arc shown in Table I. It 
is seen that the constraint on Tr is much more stringent 
than that obtained for gravitino without hadronic decay. 
Note added: While finalizing this letter, we found the 
paper by K. Jedamzik [20] which have some overlap with 
our analysis. 
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